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Abstract 

Jet photoproduction, where the two highest transverse energy (E^ 1 ) jets have E 3 ^ 1 
above 6 GeV and a jet-jet invariant mass above 23 GeV, has been studied with the ZEUS 
detector at the HERA ep collider. Resolved and direct photoproduction samples have 
been separated. The cross section as a function of the angle between the jet-jet axis 
and the beam direction in the dijet rest frame has been measured for the two samples. 
The measured angular distributions differ markedly from each other. They agree with 
the predictions of QCD calculations, where the different angular distributions reflect the 
different spins of the quark and gluon exchanged in the hard subprocess. 
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1 Introduction 



Jet photoproduction at HERA has been used to investigate various aspects of quantum chromo- 
dynamics (QCD) and the structure of the photon and the proton [1-6]. QCD calculations of jet 
cross sections can be factorized into two parts: the parton distributions in the beam particles 
and the matrix elements of the partonic hard scattering. The selection of the kinematic region 
and variables in which dijet cross sections are studied can enhance the sensitivity of the data 
either to the matrix elements or to the parton distributions. 

In a previous analysis || dijet cross sections were measured in the regime where the difference 
between the pseudorapidities^ of the two jets of highest transverse energy (E^) is small (\rf etl — 
^jei2| < 0.5). This constrained 9*, the angle between the jet-jet axis and the beam axis in 
the dijet centre of mass system, to be close to 90°. The cross section as a function of fj = 
(j]i etl -f ri^ et2 )/2 was then sensitive to the parton distributions in the photon and proton. In 
contrast, in the present analysis a cut is made on fj, resulting in a cos#* distribution which is 
sensitive to the parton dynamics. The variable cos 9* is calculated as 



Only the absolute value of cos#* can be determined because the outgoing jets are indistinguish- 
able. Measuring the distribution in cos#* is preferable to measuring the jet pseudorapidity dis- 
tribution because cos 9* is approximately invariant under the different boosts along the beam 
axis arising from the spectrum of incoming parton momenta. This minimizes the sensitivity of 
the cross section to the momentum density distribution of the partons in the beam particles. 

In leading order (LO) QCD, two types of processes lead to the photoproduction of jets. In 
direct processes (Fig. la) the photon participates in the hard scatter via either boson-gluon 
fusion or QCD Compton scattering. These processes involve a quark propagator in the s, t 
or u channel, with t and u channel processes dominating. In resolved processes (Fig. lb) the 
photon acts as a source of quarks and gluons, and only a fraction of its momentum participates 
in the hard scatter. In this case the dominant subprocesses, e.g. qg — > qg, gg — > gg and 
qq — > qq, have t-channel gluon exchange diagrams. The angular dependence of the cross section 
for resolved processes with a spin-1 gluon propagator is approximately oc (1 — | cos#*|) -2 (as 
in Rutherford scattering). This cross section rises more steeply with increasing | COS6 1 *! than 
that for direct processes with a spin-| quark propagator, where the angular dependence is 
approximately oc (1 — | cos#*|) _1 . After inclusion of all LO diagrams, QCD predicts that the 
angular distribution of the outgoing partons in resolved processes will be enhanced at high 
| cos 9* | with respect to direct photon processes. This property is expected to be preserved in 
next-to-leading order (NLO) calculations |7|]. In addition to depending upon the incoming flux 
of partons, this prediction is sensitive to the relative colour factors for each subprocess and to 
the spins of the quark and the gluon. 

The separation between direct and resolved photoproduction is only well defined at LO. To be 
able to make a measurement which can be compared to calculations at any order, the variable 

1 77 = — ln(tan |) where 9 is the polar angle with respect to the z axis, which in the ZEUS coordinate system 
is defined to be the proton direction. 




(1) 
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x obs j g usec l to define these two kinematic regions || . The variable x® BS is the fraction of the 
photon energy contributing to the production of the two highest Ej? jets. It is defined as 



= r±_ ^ (2 ) 



where E 1 is the initial photon energy. Direct processes as defined at leading order have high 
x° BS since all the photon's energy participates in the production of the hard jets, while resolved 
processes as defined at leading order have low x° BS since part of the photon's energy goes into 
the photon remnant. Thus the different spins of the quark and gluon propagators that are 
dominant in the high x® BS and the low x^ BS processes respectively should be reflected in the 
cos 8* distributions of the two samples. Henceforth throughout the following, in both the data 
and the calculations, direct and resolved samples are defined in terms of a cut on x° BS rather 
than in terms of the LO diagrams, unless explicitly stated otherwise. 

Measurements of dijet angular distributions in pp events have shown good agreement with the 
predictions of perturbative QCD |J in both fermionic and bosonic exchange processes ||. The 
measurement presented here tests QCD in a different kinematic regime and in a different pro- 
cess, and provides an opportunity to study the parton dynamics of two distinct dijet production 
mechanisms in the same experiment. 



2 Experimental Setup 

In 1994 HERA provided 820 GeV protons and 27.5 GeV positrons colliding in 153 bunches. 
Additional unpaired positron and proton bunches circulated to allow monitoring of the back- 
ground from beam-gas interactions. Events from empty bunches were used to estimate the 
background from cosmic rays. The total integrated luminosity used in this analysis is 2.6 pb _1 
collected during this running period. 



Details of the ZEUS detector have been described elsewhere |H|. The primary components used 



in this analysis are the central calorimeter and the central tracking detectors. The uranium- 
scintillator calorimeter (Tl| covers 99.7% of the total solid angle and is subdivided into three 
parts, forward (FCAL) covering 4.3 > 7] > 1.1, barrel (BCAL) covering the central region 
1.1 > T) > —0.75 and rear (RCAL) covering the backward region —0.75 > i] > —3.8 for 
an event at the nominal interaction point. Each part consists of an electromagnetic section 
followed by a hadronic section, with cell sizes of approximately 5 x 20 cm 2 (10 x 20 cm 2 in 
the rear calorimeter) and 20 x 20 cm 2 respectively. The central tracking system consists of a 
vertex detector |]T2[ and a central tracking chamber enclosed in a 1.43 T solenoidal magnetic 



field. A lead-scintillator photon calorimeter is used to measure the luminosity via the positron 
proton Bremsstrahlung process. This calorimeter is installed 100 m along the HERA tunnel 
from the interaction point in the positron direction and subtends a small angle at the interaction 



vertex fl4| . Positrons scattered through small angles are detected in a similar lead-scintillator 



calorimeter. 
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3 Event Selection 



The ZEUS detector uses a three level trigger system. At the first level events were triggered 
on a coincidence of a regional or transverse energy sum in the calorimeter and a track from 
the interaction point measured in the central tracking chamber. At the second level at least 
8 GeV total transverse energy, excluding the eight calorimeter towers immediately surrounding 
the forward beampipe, was required, and cuts on calorimeter energies and timing were used to 
suppress events caused by interactions between the proton beam and residual gas in the beam 



pipe ||15|| . At the third level, jets were found using the calorimeter cell energies and positions as 
input to a cone algorithm [[RJ. Events were required to have at least two jets of Ej? > 3.5 GeV 
and rf et < 3.0. Additional tracking cuts were made to reject proton beam-gas interactions and 
cosmic ray events. 

Further cuts are applied offline. Charged current events are rejected by a cut on the missing 
transverse momentum measured in the calorimeter. To reject remaining beam-gas and cosmic 
ray background events, tighter cuts using the final vertex position, other tracking information 
and timing information are applied. Two additional cuts are made |4j], based upon different 
measurements of the inelasticity, defined in the ZEUS frame as y = 1 — — cos8' e ) where 

E e is the energy of the incoming positron and E' e and 8' e are the energy and polar angle of the 
outgoing positron. 



Events with a positron candidate in the uranium calorimeter are removed if y e < 0.7, 
where y e is the value of y as measured assuming the positron candidate is the scattered 
positron. 



A cut is made on the Jacquet-Blondel measurement of y P7| , y,js = J2i(Ei — E zi )/2E, 



where E zi = Ei cos 9i, and Ei is the energy deposited in the calorimeter cell i which has 
a polar angle $i with respect to the measured z- vertex of the event. The sum runs over 
all calorimeter cells. For any event where the scattered positron entered the uranium 
calorimeter and either was not identified or gave y e above 0.7, the value of yjs will be 
near to one. Proton beam-gas events will have low values of y,jB- To further reduce the 
contamination from both these sources, it is required that 0.15 < yjs < 0.7. This range 
corresponds to the true y range of 0.25 < y < 0.8. 



These cuts restrict the range of the photon virtuality to less than ~ 4 GeV 2 , with a median of 
around 10~ 3 GeV 2 , which excludes deep inelastic scattering (DIS) events. 

To select dijet candidates, a cone algorithm is applied to the calorimeter cells using a cone radius 
R = J {5(f)) 2 + (Sr/) 2 = 1, where Sr](S(j)) is the distance in 7/(0) of the centre of a calorimeter 
cell from the jet centre. The jet energy measured in the ZEUS detector has been corrected 
using the Monte Carlo (MC) events described in the next section. The energy response of the 
ZEUS calorimeter is estimated by comparing jets found in the hadronic final state of the MC 
generator to jets found in the simulated detector. The average shift in jet energies is —17% and 
varies between —10% and —25% depending upon r/ je *. The description of the jet energy shift 
in the MC has been checked using tracking and the incident photon energy deduced from the 
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energy of the electron which is measured in the small-angle electron tagger. This description 
is accurate to within 5% |5j, |J. After this correction, events are required to have at least two 
jets with Ej?* > 6 GeV and rf et < 2.5. For events with three or more jets, the two highest Ej? 
jets are used to calculate all jet-related event properties. This procedure is also employed later 
in all the theoretical and MC predictions shown. 

For a given centre-of-mass energy, events at high | cos 9* | have smaller scattering angles and 
thus lower E^ 1 . In order to study the | cos#*| distribution up to | cos#*| = 0.85 without bias 
from the E^ 1 requirement, a cut of Mjj > 23 GeV has been applied. Mjj is the dijet invariant 
mass calculated using the relationship 



Mjj = V 2E J T etl E J T et2 [cosh(^' e * 1 - rp* 2 ) - cos(^' e * 1 - <^ e * 2 )] , (3) 

where je * is the azimuthal angle of the jet in the HERA frame. For two jets back-to-back in 
and with equal E J T et , Mjj « 2E J T et /^/l - \cos6*\ 2 . 

Because of the asymmetric beam energies, the two parton centre-of-mass system is typically 
boosted to the forward direction. For resolved photon processes, this effect is bigger because 
only a fraction of the photon's momentum participates in the hard scatter. MC simulations 
show that as a result, a significant fraction of jets are at rapidities (rj > 2.5) where they are not 
well measured in the ZEUS detector. To remove this bias, a cut \fj\ < 0.5 was applied, where 
fj, the average r\ of the two jets, is a measure of the boost of the dijet scattering system in the 
HERA frame: 

1 E x 

V ~ Vboost = ^H-ET- 1 ) ( 4 ) 

where x p and a; 7 are the momentum fractions of the incoming partons in the proton and photon 
respectively and E p is the incoming proton energy. In a simple (LO) 2^2 scatter, the dijet 
invariant mass is related to x p and x 1 by 



Mjj = J AE 1 E p x~ j x p = 4E e E p yx~fXp. (5) 

The requirement that the dijet system has high mass and small boost selects events with 7p 
centre-of-mass energies mostly above 190 GeV, and suppresses events with low x^ BS . Thus 
the incoming partons have approximately equal and opposite momenta and energies sufficient 
to produce dijets in a region of good acceptance for a wide range of scattering angles. In the 
range of x 7 and x p selected by these cuts the photon and proton parton distributions are fairly 
well determined by previous measurements 0, |^, ||, [L5|, |T8|| . 

After these cuts a total of 4964 events remain with I cos#*| < 0.85, of which 1982 have x° BS < 
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0.75 and 2982 have x° BS > 0.75. Events with a third jet which passes the E^ and ?7 je * cuts 
comprise 9% of the final sample. 21% of all events have their scattered positrons detected in the 
small-angle tagger, as is expected for photoproduction events ||. The beam-gas background, 
measured using unpaired bunches, is less than 1%. The contamination from events with photon 
virtualities greater than 4 GeV 2 is estimated to be 3.3% using simulated deep inelastic scattering 
(DIS) events. 
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4 Results and Discussion 



In Fig. 2 the ZEUS dijet data (black dots) are compared to the results of two QCD MC simula- 
tion programs, HERWIG58 @ (solid lines) and PYTHIA57 |Q (broken lines). The GRV [gl 



parton distributions are used for the photon and the MRSA |22[ parton distributions are used 
for the proton. The simulation programs are based on LO QCD calculations for the hard scatter 
and include parton showering and hadronization effects. The minimum transverse momentum 
of the hard scatter is set to 2.5 GeV. For both programs the direct and resolved photon pro- 
cesses are generated separately and combined according to the ratio of their generated cross 
sections. The combined MC distributions are then normalized to the number of data events. 
In the case of HERWIG resolved processes, multiparton interactions are included |2j|, as this 



has been shown to improve the simulation of the energy flow around the jet axis in low-M. 
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events |2J. All the MC events are passed through a detailed simulation of the ZEUS detector 
and the same jet energy correction procedure as was applied to the data. Fig.2a shows the 
distribution of events in Ej? , which falls steeply with increasing Ef. At low Ef the dijet 
mass cut causes the distribution to turn over. Fig. 2b shows the Mjj distribution, which falls 
steeply with increasing Mjj and extends to masses of 60 GeV. The fj distribution is shown in 
Fig. 2c, and rises with increasing fj due to the asymmetric momenta of the photon and proton. 
The fj cut leads to an t]^ et distribution which is peaked in the forward and rear directions, as 
shown in Fig. 2d. Note also that due to the fj cut, the absolute value of rf ei is restricted to be 
below 1.8. For all these distributions the data are reasonably well described by PYTHIA and 
HERWIG, although HERWIG gives more jets in the central region in r] jet and at high E 3 ^ than 
are seen in the data. 

Fig. 3a shows the uncorrected x® s distribution, where x® BS is the fraction of the proton's 
energy contributing to the production of the two highest jets: 



E Tl e^ M + E T2 e^ 



2E P 

The events lie in the range 0.006 < x° BS < 0.06, with the most probable value at Xp BS ~ 0.016. 
Fig. 3b shows the uncorrected x° BS distributionF]. Contributions from the direct photon LO 
diagrams (e.g. Fig. la) (broken line) and the resolved photon LO diagrams (e.g. Fig. lb) (dotted 
line) of HERWIG are shown together with the combined distribution (solid line). PYTHIA 
events, not shown in the figure, have a similar distribution. The data exhibit a slightly lower 
mean x° BS than the MC. A cut is applied at x° BS = 0.75 to define the direct and resolved 
photon samples. The effect of the imperfect description of the data by MC on the evaluation of 
the cross sections is evaluated by varying the x^ BS cut , and is small. The transverse energy flow 
around the jets is shown in Fig.3c for resolved and in Fig. 3d for direct events. In this sample, 
which has high Mjj values, both the HERWIG sample with multiple interactions included and 
the PYTHIA sample without multiple interactions describe the jet profiles reasonably well in 
both direct and resolved events. The requirements of high mass and small boost remove the 



2^,obs wag measurec j using the uncorrected jet transverse energies. This takes advantage of cancellations 
y G * and E-y 



between energy losses in Ef and E 1 — yjsE, 
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disagreement in the forward energy flow between data and the simulations which has been 
reported elsewhere || |], U in hard photoproduction at HERA. 

The resolution of the kinematic variables has been studied by comparing in the MC simulation 
jets reconstructed from final state particles (hadron jets) with jets reconstructed from the 
energies measured in the calorimeter (calorimeter jets), and by comparing yjg with the true y. 
The difference in | cos8*\ between the hadron and detector levels has a mean of zero, a width of 
0.03 and is approximately independent of | cos^*!, Mjj and fj. The resolution of x*? BS is 8.7% 
at x° BS = 0.75. For the variables E 3 ^ \ Mjj and y, the resolutions are 13.2%, 11.1% and 11.0% 
respectively at the values at which the cuts are applied. 

The MC samples have been used to correct for the efficiency of the trigger and selection cuts 
and migrations, and for the contamination from DIS events. The final bin-by-bin correction 
factors are approximately independent of | cos#*| and are around 1.1 for the resolved and 1.5 
for the direct cross section. These correction factors are calculated as the ratio of the purity 
(= N tTue ^ ec / N rec ) to the efficiency (= N tTuetrec /N true ) in each bin. N tTue is the number of events 
generated in the bin, N iec is the number of events reconstructed in the bin after detector 
smearing and all experimental cuts, and N trne ^ Tec is the number of events which were both 
generated in the bin and reconstructed in that bin. The difference between the correction 
factors for resolved and direct events is due to the fact that the E^ 1 and Mjj distributions are 
steeper in resolved processes than in the direct, and thus migrations from below the cuts are 
more significant, giving a lower purity for resolved than direct. 

The sensitivity of the measured cross sections to the selection cuts has been investigated by 
varying the cuts on the reconstructed variables in the data and HERWIG MC samples and 
re-evaluating the cross section. In addition, the cross section was re-evaluated using a different 
ratio of the direct and resolved contributions, and using the PYTHIA sample. The systematic 
uncertainty arising from the subtraction of the contamination coming from DIS was estimated 
by using two different DIS Monte Carlo generators and two different positron-finding algo- 



rithms [Tq, 18 1 . We have also allowed for the possibility that the detector simulation may 



overestimate the detector energy response by up to 5%, as mentioned in section 3. 

The cross sections da/d\ cos#*| for ep — > dijets +X in the range Mjj > 23 GeV, E 3 ^ 1 > 6 GeV, 
\fj\ < 0.5, 0.25 < y < 0.8 and for virtualities of the exchanged photon < 4 GeV 2 are listed in 
table 1. For the direct photoproduction sample the ep cross section rises from around 0.8 nb 
at | cos#*| = to around 6 nb at | cos#*| = 0.85. For the resolved photoproduction sample the 
ep cross section rises from around 0.2 nb at | cos#*| = to around 4 nb at | cos#*| = 0.85. In 
Fig. 4a and b the two cross sections have been normalized to unity at | cos 6* \ = so that the 
shapes can be compared directly. To determine the normalization, the 5 lowest | cos 6* | bins 
were fitted to the functional form expected for the dominant cross section in each sample (as 
described in section 1), which is slowly rising in this region. The only free parameter in this 
fit is the normalisation of the function. The intercept of the function then gives the factor by 
which the data are divided. In Fig. 4a the data are compared to LO and partial [] NLO QCD 
parton level calculations @. The resolved cross section is seen to rise more steeply than the 
direct cross section with increasing | cos 8* \ . There is good agreement between data and theory, 

3 The NLO calculations do not yet include terms of 0(a 3 x /), where / is the flux of partons in the photon. 
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verifying the expected effects of the spins of the quark and gluon propagators. This is also seen 
in Fig. 4b, where the data are compared to the HERWIG and PYTHIA predictions. For both 
of the generators, MRSA(proton) and GRVLO (photon) |2]J parton parameterisation sets are 
used. The MC curves were also made with different parton distribution sets in the photon and 
proton. The largest change seen was for the LAC1 J^J photon parton distribution set, where 
there is a 35% increase in the prediction in the highest | cos#*| bin. In all the calculations, 
resolved and direct samples are defined on the basis of the cut on x® BS , exactly as in the data. 
Both the MC and the parton level calculations show the same behaviour, suggesting that the 
effects of parton showers and hadronisation on this distribution are small. The data are in good 
agreement with the MC and with the calculations. 

The only difference in the cuts applied to the two samples is the cut on x° BS which is used 
to separate the resolved and direct samples. A MC study has been performed to investigate 
whether the different behaviour of the two cross sections might be a kinematic effect of this 
cut. The study shows that the tail of the resolved photon LO diagrams (e.g. Fig. lb) which lie 
above x° BS = 0.75 retain the steeply rising | cos^*! distribution characteristic of the low-x® BS 
events. Thus the different behaviour of the two cross sections is due to the different partonic 
subprocesses which contribute in the different x° BS regions. 



5 Summary and Conclusion 

The dijet angular cross section da/d\ cos 9* | , where 9* is the jet scattering angle in the dijet c.m.s, 
has been measured for ep — > 2 or more jets with transverse energies E^ 1 > 6 GeV, dijet invariant 
mass Mjj > 23 GeV and average pseudorapidity \fj\ < 0.5, in the range 0.25 < y < 0.8 and for 
virtualities of the exchanged photon less than 4 GeV 2 . The cross section has been measured 
for resolved (x° BS < 0.75) and direct (x° BS > 0.75) processes. The angular dependence for 
the two samples is significantly different. 

The dependence of da/d\ cos#*| on | cos^*! reflects the different spins of the quark and gluon 
propagators and the relative contributions of the underlying subprocesses in resolved and direct 
photoproduction. In LO QCD the cross section rises faster with increasing | cos#*| for resolved 
photoproduction, where processes involving spin-1 gluon exchange dominate, than for direct 
photoproduction, where processes involving spin-1/2 quark exchange dominate. These expec- 
tations are preserved in NLO QCD calculations and in Monte Carlo simulations which include 
parton showering and hadronisation models. The | cos 9* | dependence of the measured cross 
sections is in good agreement with these theoretical predictions and thus confirms fundamental 
aspects of quantum chromo dynamics. 
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COS0*\ 






Direct (nb / unit cos 








range 


da/d cos 9* 


± stat. 


syst. error 


PYTHIA 


E scale -5% 


E scale ±5% 


0.000 to 0.085 


0.80 ±0.12 


±0.06/-0.09 


-0.08 




-0.07 


0.06 


0.085 to 0.170 


0.83 ±0.13 


+0.07/-0.13 


-0.10 




-0.13 


0.11 


0.170 to 0.255 


0.98 ±0.15 


+0.02/-0.13 


-0.17 




-0.16 


0.16 


0.255 to 0.340 


0.90 ±0.14 


±0.02/-0.08 


-0.13 




-0.12 


0.12 


0.340 to 0.425 


1.26 ±0.19 


+0.15/-0.16 


-0.30 




-0.21 


0.26 


0.425 to 0.510 


1.50 ±0.18 


+0.08/-0.19 


-0.26 




-0.26 


0.24 


0.510 to 0.595 


1.41 ±0.17 


+0.18/-0.08 


-0.10 




-0.21 


0.14 


0.595 to 0.680 


2.70 ±0.28 


±0.36/-0.35 


-0.27 




-0.53 


0.51 


0.680 to 0.765 


3.81 ±0.33 


+0.14/-0.29 


-0.30 




-0.57 


0.36 


0.765 to 0.850 


5.99 ±0.48 


+0.81/-0.40 


-0.37 




-0.96 


0.86 


| cos 0*| 


Resolved (nb / unit cos#*|) 


range 


da/d cos 6* 


± stat. 


syst. error 


PYTHIA 


E scale -5% 


E scale ±5% 


0.000 to 0.085 


0.21 ±0.08 


±0.06/-0.08 


-0.09 




-0.04 


0.00 


0.085 to 0.170 


0.23 ±0.12 


+0.09/-0.12 


-0.08 




-0.06 


-0.02 


0.170 to 0.255 


0.23 ±0.10 


±0.04/-0.07 


0.01 




-0.06 


0.12 


0.255 to 0.340 


0.28 ±0.09 


+0.13/-0.02 


0.00 




-0.02 


0.04 


0.340 to 0.425 


0.27 ±0.10 


+0.10/-0.06 


-0.03 




0.01 


0.02 


0.425 to 0.510 


0.42 ±0.11 


+0.11/-0.04 


-0.03 




-0.06 


0.11 


0.510 to 0.595 


0.77 ±0.17 


±0.08/-0.08 


0.21 




-0.13 


0.10 


0.595 to 0.680 


1.13 ±0.21 


+0.16/-0.26 


-0.02 




-0.22 


0.07 


0.680 to 0.765 


2.03 ±0.25 


+0.36/-0.19 


0.01 




-0.39 


0.63 


0.765 to 0.850 


4.31 ±0.43 


±0.47/-0.39 


0.15 




-0.85 


0.97 



Table 1: Differential cross sections da/d \ cos 0*\ for ep — > 2 or more jets with E$T > 6 GeV, 
Mjj > 23 GeV and \fj\ < 0.5, in the range 0.25 < y < 0.8 and for virtualities of the exchanged 
photon less than 4 GeV 2 , for direct and resolved processes. The first column shows the cross 
section and statistical errors in units of nb, the second column the uncorrelated systematic 
errors, the third column the shift in the cross section when the correction is evaluated using 
PYTHIA instead of HERWIG and the fourth and fifth columns the shifts in the cross section 
when the calorimeter energy calibration is varied by ±5%. There is an additional overall 
uncertainty of 3.3% arising from the luminosity measurement. 
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ZEUS 1994 




Figure 2: (a) E^ 1 , (b) Mjj, (c) fj and (d) if et distributions. Raw ZEUS data are compared to 
the results of the HERWIG58 (solid line) and PYTHIA57 (broken line) Monte Carlo models 
after simulation of all detector effects and application of selection cuts. The MC samples have 
been normalized to the number of events in the data. Only statistical errors are shown. 
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Figure 3: Raw ZEUS data compared to results of Monte Carlo models after simulation of 
detector effects and application of selection cuts, (a) x® BS and (b) x^ B distribution. The MC 
samples have been normalized to the number of events in the data. Figures (c) and (d) show 
the uncorrected transverse energy flow l/NdE T /dSr) around the jet axis, for cells within one 
radian in of the jet axis, for (c) resolved (x° BS < 0.75), (d) direct (x° BS > 0.75) events. Only 
statistical errors are shown. In (a), (c) and (d), the histograms are HERWIG58 (solid line) and 
PYTHIA57 (broken line) In (b), the solid line is HERWIG58, the broken line is the distribution 
of the LO direct diagrams (e.g. Fig. la) and the dotted line is that of the LO resolved diagrams 
(e.g. Fig.lb). 
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Figure 4: da/d\ cos#*| normalized to one at cos 9* = for resolved (black dots) and direct (open 
circles) photoproduction. In (a), the ZEUS data are compared to the NLO prediction (solid 
line) and the LO prediction (broken line). The parton distribution sets used in the calculation 
are CTEQ3M [26] for the proton and GRV (LO) [21] for the photon. In (b), the broken line is 
the PYTHIA distribution and solid line is HERWIG distribution. The inner error bars are the 
statistical errors, the outer error bars are the sum in quadrature of the statistical and systematic 
errors, excluding the energy scale and luminosity uncertainties (see Table 1). 
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